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Abstract 

We study the nonresonant part of the ''Be(p, 7)®B reaction using a three- 
cluster resonating group model that is variationally converged and virtually 
complete in the ^He+^He+p model space. The importance of using adequate 
nucleon-nucleon interaction is demonstrated. We find that the low-energy 
astrophysical S-factor is linearly correlated with the quadrupole moment of 
^Be. A range of parameters is found where the most important ''Be and 
''Li properties are reproduced simultaneously; the corresponding S'-factor at 
^cm = 20 keV is 24.6 - 26.1 eV-b. 

PACS numbers: 25.40.Lw, 21.60.Gx, 27.20.-hn, 95.30.Cq 
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The flux of high-energy neutrinos generated in the solar core is directly proportional to 
the "^Be^j, 7)^B reaction rate. Thus, knowledge of Sn, the ''Be(p, 7)^8 5'-factor at solar 
energies (center-of-mass energy 20 keV), is crucial to conclusions drawn from present 
(Homestake, Kamiokande) and future (SNO, Superkamiokande) solar neutrino experiments 
Despite extensive experimental efforts, the ^Be(p, 7)^8 cross section is still the most 
uncertain nuclear input to the standard solar model |jl|J^, due to a signiflcant spread among 
the values of 5*17 deduced from the various experiments (direct capture Sn = 18 — 28 
eV-b and Coulomb break-up 0: 6*17 = 16.7 ± 3.5 eV-b). Theoretical estimates also vary 
(5*17 = 16 — 30 eV-b) making these predictions rather unreliable. 

The aim of this Letter is to constrain more tightly the theoretical value of 5*17. To 
this end, we study the ^Be(p, 7)^8 reaction in a microscopic three-cluster (^He+^He-|-p) 
approach. This model is currently the closest approximation to a full solution of the micro- 
scopic eight-nucleon problem with a consistent treatment of bound and scattering states. As 
we will demonstrate below, our approach is superior to all previous studies of the low-energy 
^Be(p, 7)^B reaction. 

Adopting a microscopic three-cluster (^He-|-^He-|-p) ansatz for the eight-nucleon system, 
our trial function reads 



{ij)k,S,li,h,L 

where the indices i,j, and k denote any of the labels ^He, ^He, and p. In (||) A is the inter- 
cluster antisymmetrizer, the cluster internal states $ are translationly invariant harmonic 
oscillator shell model states, the p vectors are the intercluster Jacobi coordinates, and [. . .] 
denotes angular momentum coupling. In the sum over S,li,l2, and L we include all angu- 
lar momentum conflgurations of any signiflcance. This same model was used in in the 
study of the ground state of ^B; further details on the model space and other aspects can be 
found there. The intercluster dynamics is determined by inserting (|I]) into the eight-nucleon 
Schrodinger equation using the two-nucleon strong and Coulomb interactions. In addition to 
the full model space calculation, which contains all three possible arrangements of the three 
clusters, we also present a restricted calculation involving only (^He,^He)p conflgurations 
(^Be-|-p type model space), analogous to simple ^Be+p potential model studies, e.g. Pj. 

It is well known that the low-energy ^Be(p, 7)^B cross section is strongly dominated by 
El capture. (Previous microscopic calculations have shown that Ml capture only plays a 
role in the vicinity of the l"*" resonance at E = 640 keV and is negligible at astrophysical 
energies 0, while E2 capture is tiny at < 500 keV and can safely be ignored.) We have 
therefore calculated the El capture cross section into the *B ground state in perturbation 
theory (as outlined for example in |^), describing the initial scattering states and the ^B 
ground state by the many-body wave functions determined in our microscopic three-cluster 
approach. 

The capture cross section contains the bound '^B and the scattering ''Be-l-p wave func- 
tions. At low energies, deep below the Coulomb barrier, the capture takes place at large 
^Be-p distances, which means that these wave functions must be accurate to distances of 
a few hundred fermis, which requires a reliable method to determine the unknown relative 
motion functions x in (0). We expand these functions in terms of products of basis functions 
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of the Jacobi coordinates, which allow us to reduce the three-cluster wave functions (|l|) to 
equivalent two-cluster forms For example, if we have N basis functions between the 



"^He and ^He clusters in the ('^He,'^He)p partition, after this reformulation we get a ^Be+p 
two-cluster wave function with the ^Be ground state and N — 1 (unphysical, continuum) 
excited states. 

We use the variational Siegert method to determine the ®B boundstate [|lT|. The trial 



state contains tempered Gaussian functions [0 plus a term with the correct outgoing Whit- 
taker assymptotics, W^^\kp), in the "^Be+p partitions. (Here k is the wave number cor- 
responding to the *B binding energy relative to the ^Be+p threshold.) Using such a trial 
function in a linear variational method leads to a transcendental equation for the binding 
energy, which can be solved iteratively. To be able to calculate every many-body matrix 
element analytically, we match the external Whittaker functions with internal Gaussians, 



using a modified version of the technique described in |13|. The accuracy of this procedure is 



better than 1-2% in 5*17. For comparison, the uncertainity coming from the different possible 
ways of handling the proton-neutron mass difference is estimated to be 2.5% [0]. 

The scattering wave functions were calculated using the variational Kohn-Hulthen 



method |]T3[ , which ensures the correct scattering assymptotics. To achieve high accuracy we 
avoid the use of complex wave functions and so neglect channel coupling between different 
angular momentum channels. Note that this approximation is certainly justified at astro- 
physical energies where the capture occurs far outside the range of the strong forces. The 
present scattering solution is numerically well-conditioned for E > 3 keV, and its accuracy 
is better than 0.1 %. The technical details and further physical implications of the model 
will be published elsewhere |jl5|. 

The bulk of our calculations use the Minnesota (MN) effective nucleon-nucleon inter- 
action |T^, which contains central and spin-orbit terms. This force reproduces the most 
important properties of the low-energy N + N and ^He+A^ scattering phase shifts and, as 
we show below, appears to be best suited for the problem at hand. However, we will also 
compare these results to those obtained with other effective NN interactions. Note that the 
tensor component of the effective NN interaction in microscopic cluster models is not well 
constrained |^ and is usually ignored. Nevertheless we have also performed a calculation 
including a tensor force, which, at the least, gives the correct low-energy order of the triplet- 
odd N + N phase shifts 0. As is customary, we will present our results in terms of the 
astrophysical ^-factor 

S{E) = a{E)Eexp [27rr]{E)], (2) 

where 77 is the Sommerfeld parameter. 

The free parameters in our model are the size parameter [f3) in the ^He and ^He cluster 
model functions (technical reasons force us to use the same value for both ^He and ^He), 
the exchange mixture parameter of the central part of the effective NN interaction, and 
the strength of the spin-orbit force. It is generally preferable to adjust these parameters to 
independent data. However, a meaningful study of the '^Be{p,'j)^B reaction at low energies 
requires the exact reproduction of the experimental ^B binding energy (137 keV), which we 
have guaranteed by the appropriate choice of the exchange mixture parameter. The strength 
of the spin-orbit force was adjusted to the experimental splitting between the 3/2~ and l/2~ 
^Be states. We have varied j3, thus changing our description of the ^Be properties. 
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As is demonstrated by the open circles in Fig. ^ Sn scales linearly with the quadrupole 
moment of ^Be, Qtbe- This linear dependence can be understood as follows. As the capture 
process takes place at very large ^Be-p distances, where the bound state wave function is 
proportional to the fixed Whittaker function, '^a^{p) ■ p = cWjj^\kp), the low-energy cross 
section depends almost exclusively on the square of the asymptotic normalization factor 
c. Let us compare calculations with different ''Be wave functions, which give different ^Be 
radius, quadrupole moment, etc., but with fixed binding energy of ^B. The effective local po- 
tentials between ^Be and p have different radii, which means that the height of the Coulomb 
barrier is larger if the potential radius (and the ^Be radius) is smaller. Consequently, the 
probability of finding the proton in the outside region decreases as the size of the ^Be nucleus 
becomes smaller. But as the shape of the external wave function is fixed, this smaller prob- 
ability must stem from a smaller normalization constant c. It is easy to see that this leads 
(f, and consequently Su, to be linearly proportional to either r?gg or Qtbg- Note that this 
relation is not changed if a tensor component is added to the MN interaction (see triangle in 
Fig. m). We find the same linear Sn-Qr^^ relation in our truncated calculation considering 
only the ''Be+p model space. Results of these restricted calculations are shown in Fig. |I| as 
full circles. 

Unfortunately the linear relation is not sufficient to determine 5*17 indirectly by measuring 
the ''Be quadrupole moment, as this relation depends upon the effective NN interaction used. 
To demonstrate this we have performed calculations within the ^Be+p model space using the 
Volkov force V2 and the modified Hasegawa-Nagata (MHN) force, both of which have been 
used in previous microscopic cluster calculations of the ''Be(p, 7)*B reaction at low energies 
P JT7|JI8|] . While both forces also show the linear dependence between Sn and the ^Be 
quadrupole moment, the V2 force yields larger values for Su for a given Qy^^ (diamonds in 
Fig. |I|), while the MHN force yields smaller values (squares). These differences can be traced 
to the different quality of the description of the N + N systems (phase shifts, energy and 
radius of the deuteron) by these forces. For example, while the MN force well reproduces the 
experimental deuteron energy and radius, the V2 force underbinds the deuteron by 1.6 MeV 
(however, it unphysically binds the singlet dinucleon states) and the MHN force overbinds 
it by 4.4 MeV. We note that the M3Y interaction, which was used in Ref. [|19| in an external 
capture approach to predict a very small ''Be(p,7)'^B cross section {Sn =16.5 eV-b), also 
overbinds the deuteron. From this discussion we conclude that the Minnesota force is by far 
the most carefully constructed force available in the cluster literature; we will adopt it in the 
following. Relatedly, cluster calculations using the MHN and V2 forces should be regarded 
with care. 

Accepting the MN force as adequate for the eight-nucleon problem, our result for Sn 
could be read off Fig. |I| if the ^Be quadrupole moment were known. Absent this information, 
we will estimate a best 5*17 value by constraining the ^He and ^He cluster size parameter to 
reproduce (i) the binding energy of ^Be with respect to "^He+^He; (ii) the squared sum of 
the ^He and ^He radii; (iii) the quadrupole moment of ^Li (as a surrogate for the unknown 
quadrupole moment of the analog nucleus ''Be). These requirements ensure that both the ^Be 
bound states and the '^He-^He relative motion are well described. The second requirement 
is fulfilled by choosing P = 0.4 fm~^. With this choice, the ''Be ground state is slightly 
underbound by 200 keV, while the quadrupole moment of ^Li is calculated as Q = —4.10 
e-fm^, to be compared with the experimental value — 4.05±0.08 e-fm^ PO. 
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We conclude that the three requirements above can be reasonably fulfilled simultaneously. 
The corresponding Sn value is then 26.1 eV-b, while the ^Be quadrupole moment is —6.9 
e-fm^. Our approach then calculates the quadrupole moment of *B as 7.45 e-fm^, while the 
experimental value is (6.83±0.21) e-fm^ Even if one concludes from a comparison of 
our ^Li and quadrupole moments with experiment that our ''Be quadrupole moment is 
also slightly too large, we note that a 10% reduction in this quantity would only decrease 
Sn to 23.5 eV-b. 

If we use the same cluster size parameter in the restricted ^Be+p space as in the full 
calculation (/3 = 0.4 fm^), we find that the ^Be nucleus is overbound (by 600 keV), while 
its quadrupole moment is reduced to Q = —6.0 e-fm^. The quadrupole moments of ^Li 
{Q = —3.46 e-fm^) and ®B {Q = 6.55 e-fm^) are slightly smaller than the experimental 
values. In this restricted calculation we find 5*17 to be 24.6 eV-b. 

If we consider that both the full and restricted ^Be+p model spaces predict the same 
linear dependence of Sn on the ^Be quadrupole moment and that these calculations bracket 
the experimental ^Li and ^B quadrupole moments, we conclude that the microscopic three- 
cluster calculations predict 5*17 to be in the range between 24.6 and 26.1 eV-b. This does 
not support speculations that 5*17 might be noticeably smaller {Sn =16.5, 16.9, and 17 eV-b 
in |T^, ||22|, and respectively) than currently accepted in the standard solar model. We 



note, however, that 5*17 deduced from our model is consistent with the value deduced from 
the direct capture data (22±2 eV-b @, and 24±2 eV-b @], respectively). 

Less elaborate microscopic cluster calculation have been presented in Refs. P, p4| , p!8| , p!7 



While the two earlier studies P,p4| were restricted to a simple ''Be+p model space, Ref. |T7 
recently improved these studies by including a ^Li+^He rearrangement channel. However, 
in 1 17] the ^Be nucleus is described by only one Gaussian basis function between ^He and 



^He, which means that the three-cluster wave function is not free for the variational method. 
Letting the trial function more fiexible would result in the collapse of the artificially fixed 
wave function. Moreover, in []T^ the description of the ^Be nucleus is rather unphysical, 
as it is unbound relative to the ^He-|-^He threshold. In [0 there are two basis functions 
for ^Be, with carefully chosen parameters, and the most important angular momentum 
configurations of the ''Be+p type partition are present. In the present model we use six 
states for ^Be (and ten in the ^Be-p relative motion, and six in all other relative motions) 
and include all relevant angular momentum channel. Our test calculations showed that 
the present three-cluster model space is virtually complete (see also the discussion in 0), 
which means that our results do not contain the side-effects of an unconverged or incomplete 
model. Although the incompleteness of the previous works makes the comparison difficult, 
our results are qualitatively in good agreement with [|T^ and |TB]. 

In Fig. we show the energy-dependence of the S-factor, calculated with the ^Be-|-p 
model space, the MN force, and j3 = OA fm^. At low energies our calculated S-factor is in 
rather close agreement with the direct capture data of Ref. |2^. Although our ^-factor also 
agrees well with the data of for E > 1 MeV, its energy dependence might change, if 
the coupling between the different angular momentum channels in the ''Be+p scattering is 
taken into account. 

In summary, we have studied the ^Be(p, 7)®B reaction in a microscopic model. At low 
energies this model is virtually complete in the three-cluster model space. We found that 
the low-energy astrophysical S'-factor is strongly correlated with the properties of ^Be (e.g., 
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its quadrupole moment). For a set of parameters that reproduce simultaneously the most 
important properties of ''Be, ^Li and we predict Qrse to be between —6.0 e-fm^ and —6.9 
e-fm^ and find Sn = 24.6 — 26.1 eV-b, in good agreement with direct capture results and 
the currently accepted value in the standard solar model. Our calculation, thus does not 
support the recently suggested smaller values for Su If it turns out that the S- 

factor is considerably lower than our present value, then the present three-cluster approach 
is inappropriate and physics beyond our model (larger eight-body model space, improved 
effective interaction) has to be invoked. We have also shown that the NN interaction used in 
cluster models must be carefully chosen. Although we found that only the Minnesota force 
was suitable for the present work, the construction and use of other high quality interactions 
would be indespensable. We also note that a precise measurement of the ''Be quadrupole 
moment or radius could test the self-consistency of our conclusions. 

This work was supported by the Fulbright Foundation (A. C.) and NSF Grant Nos. 
PHY90-13248 and PHY91-15574 (USA), and by OTKA Grant Nos. 3010 and F4348 (Hun- 
gary). We thank Dr. E. Kolbe for useful discussions and Dr. Z. Papp for his generous help 
in numerical methods. 
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FIGURES 



FIG. 1. The astrophysical 5-factor of the ^Be(p, 7)^6 reaction as a function of the negative of 
the ''Be quadrupole moment. The symbols are explained in the text. 



FIG. 2. Energy dependence of the ^Be(p, 7)^B astrophysical 5-factor. The symbols denote the 
experimental data of Ref. (open circles), Ref. |^] (filled circles), and Ref. |27] (squares). The 
inset shows the low-energy part on a magnified scale. See the error bars in the original references. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9408001vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9408001vl 



